BATCH-GE : batch analysis of next-generation sequencing data for genome editing assessment

Targeted mutagenesis by the CRISPR/Cas9 system is currently revolutionizing genetics. The ease of this technique has enabled genome engineering in-vitro and in a range of model organisms and has pushed experimental dimensions to unprecedented proportions. Due to its tremendous progress in terms of speed, read length, throughput and cost, Next-Generation Sequencing (NGS) has been increasingly used for the analysis of CRISPR/Cas9 genome editing experiments. However, the current tools for genome editing assessment lack flexibility and fall short in the analysis of large amounts of NGS data. Therefore, we designed BATCH-GE, an easy-to-use bioinformatics tool for batch analysis of NGS-generated genome editing data, available from https://github.com/WouterSteyaert/BATCH-GE.git. BATCH-GE detects and reports indel mutations and other precise genome editing events and calculates the corresponding mutagenesis efficiencies for a large number of samples in parallel. Furthermore, this new tool provides flexibility by allowing the user to adapt a number of input variables. The performance of BATCH-GE was evaluated in two genome editing experiments, aiming to generate knock-out and knock-in zebrafish mutants. This tool will not only contribute to the evaluation of CRISPR/Cas9-based experiments, but will be of use in any genome editing experiment and has the ability to analyze data from every organism with a sequenced genome

Recessive osteogenesis imperfecta caused by LEPRE1 mutations: clinical documentation and identification of the splice form responsible for prolyl 3-hydroxylation

Abstract: Background: Recessive forms of osteogenesis imperfecta (OI) may be caused by mutations in LEPRE1, encoding prolyl 3-hydroxylase-1 (P3H1) or in CRTAP, encoding cartilage associated protein. These proteins constitute together with cyclophilin B (CyPB) the prolyl 3-hydroxylation complex that hydroxylates the Pro986 residue in both the type I and type II collagen alpha 1-chains. Methods: We screened LEPRE1, CRTAP and PPIB (encoding CyPB) in a European/Middle Eastern cohort of 20 lethal/severe OI patients without a type I collagen mutation. Results: Four novel homozygous and compound heterozygous mutations were identified in LEPRE1 in four probands. Two probands survived the neonatal period, including one patient who is the eldest reported patient (17(7/12) years) so far with P3H1 deficiency. At birth, clinical and radiologic features were hardly distinguishable from those in patients with autosomal dominant (AD) severe/lethal OI. Follow-up data reveal that the longer lived patients develop a severe osteochondrodysplasia that overlaps with, but has some distinctive features from, AD OI. A new splice site mutation was identified in two of the four probands, affecting only one of three LEPRE1 mRNA splice forms, detected in this study. The affected splice form encodes a 736 amino acid (AA) protein with a "KDEL'' endoplasmic reticulum retention signal. While western blotting and immunocytochemical analysis of fibroblast cultures revealed absence of this P3H1 protein, mass spectrometry and SDS-urea-PAGE data showed severe reduction of alpha 1(I) Pro986 3-hydroxylation and overmodification of type I (pro) collagen chains in skin fibroblasts of the patients. Conclusion: These findings suggest that the 3-hydroxylation function of P3H1 is restricted to the 736AA splice form

CRISPR/Cas9 mediated knockout of rb1 and rbl1 leads to rapid and penetrant retinoblastoma development in Xenopus tropicalis

Retinoblastoma is a pediatric eye tumor in which bi-allelic inactivation of the Retinoblastoma 1 (RB1) gene is the initiating genetic lesion. Although recently curative rates of retinoblastoma have increased, there are at this time no molecular targeted therapies available. This is, in part, due to the lack of highly penetrant and rapid retinoblastoma animal models that facilitate rapid identification of targets that allow therapeutic intervention. Different mouse models are available, all based on genetic deactivation of both Rb1 and Retinoblastoma-like 1 (Rbl1), and each showing different kinetics of retinoblastoma development. Here, we show by CRISPR/Cas9 techniques that similar to the mouse, neither rb1 nor rbl1 single mosaic mutant Xenopus tropicalis develop tumors, whereas rb1/rbl1 double mosaic mutant tadpoles rapidly develop retinoblastoma. Moreover, occasionally presence of pinealoblastoma (trilateral retinoblastoma) was detected. We thus present the first CRISPR/Cas9 mediated cancer model in Xenopus tropicalis and the first genuine genetic non-mammalian retinoblastoma model. The rapid kinetics of our model paves the way for use as a pre-clinical model. Additionally, this retinoblastoma model provides unique possibilities for fast elucidation of novel drug targets by triple multiplex CRISPR/Cas9 gRNA injections (rb1 + rbl1 + modifier gene) in order to address the clinically unmet need of targeted retinoblastoma therapy

CRISPR/Cas9-mediated homology-directed repair by ssODNs in zebrafish induces complex mutational patterns resulting from genomic integration of repair-template fragments

Targeted genome editing by CRISPR/Cas9 is extremely well fitted to generate gene disruptions, although precise sequence replacement by CRISPR/Cas9-mediated homology-directed repair (HDR) suffers from low efficiency, impeding its use for high-throughput knock-in disease modeling. In this study, we used next-generation sequencing (NGS) analysis to determine the efficiency and reliability of CRISPR/Cas9-mediated HDR using several types of single-stranded oligodeoxynucleotide (ssODN) repair templates for the introduction of disease-relevant point mutations in the zebrafish genome. Our results suggest that HDR rates are strongly determined by repair-template composition, with the most influential factor being homology-arm length. However, we found that repair using ssODNs does not only lead to precise sequence replacement but also induces integration of repair-template fragments at the Cas9 cut site. We observed that error-free repair occurs at a relatively constant rate of 1-4% when using different repair templates, which was sufficient for transmission of point mutations to the F1 generation. On the other hand, erroneous repair mainly accounts for the variability in repair rate between the different repair templates. To further improve error-free HDR rates, elucidating the mechanism behind this erroneous repair is essential. We show that the error-prone nature of ssODN-mediated repair, believed to act via synthesis-dependent strand annealing (SDSA), is most likely due to DNA synthesis errors. In conclusion, caution is warranted when using ssODNs for the generation of knock-in models or for therapeutic applications. We recommend the application of in-depth NGS analysis to examine both the efficiency and error-free nature of HDR events

Slc2a10 knock-out mice deficient in ascorbic acid synthesis recapitulate aspects of arterial tortuosity syndrome and display mitochondrial respiration defects

Arterial tortuosity syndrome (ATS) is a recessively inherited connective tissue disorder, mainly characterized by tortuosity and aneurysm formation of the major arteries. ATS is caused by loss-of-function mutations in SLC2A10, encoding the facilitative glucose transporter GLUT10. Former studies implicated GLUT10 in the transport of dehydroascorbic acid, the oxidized form of ascorbic acid (AA). Mouse models carrying homozygous Slc2a10 missense mutations did not recapitulate the human phenotype. Since mice, in contrast to humans, are able to intracellularly synthesize AA, we generated a novel ATS mouse model, deficient for Slc2a10 as well as Gulo, which encodes for L-gulonolactone oxidase, an enzyme catalyzing the final step in AA biosynthesis in mouse. Gulo;Slc2a10 double knock-out mice showed mild phenotypic anomalies, which were absent in single knock-out controls. While Gulo;Slc2a10 double knock-out mice did not fully phenocopy human ATS, histological and immunocytochemical analysis revealed compromised extracellular matrix formation. Transforming growth factor beta signaling remained unaltered, while mitochondrial function was compromised in smooth muscle cells derived from Gulo;Slc2a10 double knock-out mice. Altogether, our data add evidence that ATS is an ascorbate compartmentalization disorder, but additional factors underlying the observed phenotype in humans remain to be determined

Zebrafish : a resourceful vertebrate model to investigate skeletal disorders

Animal models are essential tools for addressing fundamental scientific questions about skeletal diseases and for the development of new therapeutic approaches. Traditionally, mice have been the most common model organism in biomedical research, but their use is hampered by several limitations including complex generation, demanding investigation of early developmental stages, regulatory restrictions on breeding, and high maintenance cost. The zebrafish has been used as an efficient alternative vertebrate model for the study of human skeletal diseases, thanks to its easy genetic manipulation, high fecundity, external fertilization, transparency of rapidly developing embryos, and low maintenance cost. Furthermore, zebrafish share similar skeletal cells and ossification types with mammals. In the last decades, the use of both forward and new reverse genetics techniques has resulted in the generation of many mutant lines carrying skeletal phenotypes associated with human diseases. In addition, transgenic lines expressing fluorescent proteins under bone cell- or pathway- specific promoters enable in vivo imaging of differentiation and signaling at the cellular level. Despite the small size of the zebrafish, many traditional techniques for skeletal phenotyping, such as x-ray and microCT imaging and histological approaches, can be applied using the appropriate equipment and custom protocols. The ability of adult zebrafish to remodel skeletal tissues can be exploited as a unique tool to investigate bone formation and repair. Finally, the permeability of embryos to chemicals dissolved in water, together with the availability of large numbers of small-sized animals makes zebrafish a perfect model for high-throughput bone anabolic drug screening. This review aims to discuss the techniques that make zebrafish a powerful model to investigate the molecular and physiological basis of skeletal disorders

Purification of high-quality RNA from a small number of fluorescence activated cell sorted zebrafish cells for RNA sequencing purposes

Background: Transgenic zebrafish lines with the expression of a fluorescent reporter under the control of a cell-type specific promoter, enable transcriptome analysis of FACS sorted cell populations. RNA quality and yield are key determinant factors for accurate expression profiling. Limited cell number and FACS induced cellular stress make RNA isolation of sorted zebrafish cells a delicate process. We aimed to optimize a workflow to extract sufficient amounts of high-quality RNA from a limited number of FACS sorted cells from Tg(fli1a:GFP) zebrafish embryos, which can be used for accurate gene expression analysis. Results: We evaluated two suitable RNA isolation kits (theRNAqueous micro and the RNeasy plus micro kit) and determined that sorting cells directly into lysis buffer is a critical step for success. For low cell numbers, this ensures direct cell lysis, protects RNA from degradation and results in a higher RNA quality and yield. We showed that this works well up to 0.5x dilution of the lysis buffer with sorted cells. In our sort settings, this corresponded to 30,000 and 75,000 cells for the RNAqueous micro kit and RNeasy plus micro kit respectively. Sorting more cells dilutes the lysis buffer too much and requires the use of a collection buffer. We also demonstrated that an additional genomic DNA removal step after RNA isolation is required to completely clear the RNA from any contaminating genomic DNA. For cDNA synthesis and library preparation, we combined SmartSeq v4 full length cDNA library amplification, Nextera XT tagmentation and sample barcoding. Using this workflow, we were able to generate highly reproducible RNA sequencing results. Conclusions: The presented optimized workflow enables to generate high quality RNA and allows accurate transcriptome profiling of small populations of sorted zebrafish cells

Maximizing CRISPR/Cas9 phenotype penetrance applying predictive modeling of editing outcomes in Xenopus and zebrafish embryos

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Naert, T., Tulkens, D., Edwards, N. A., Carron, M., Shaidani, N. I., Wlizla, M., Boel, A., Demuynck, S., Horb, M. E., Coucke, P., Willaert, A., Zorn, A. M., & Vleminckx, K. Maximizing CRISPR/Cas9 phenotype penetrance applying predictive modeling of editing outcomes in Xenopus and zebrafish embryos. Scientific Reports, 10(1), (2020): 14662, doi:10.1038/s41598-020-71412-0.CRISPR/Cas9 genome editing has revolutionized functional genomics in vertebrates. However, CRISPR/Cas9 edited F0 animals too often demonstrate variable phenotypic penetrance due to the mosaic nature of editing outcomes after double strand break (DSB) repair. Even with high efficiency levels of genome editing, phenotypes may be obscured by proportional presence of in-frame mutations that still produce functional protein. Recently, studies in cell culture systems have shown that the nature of CRISPR/Cas9-mediated mutations can be dependent on local sequence context and can be predicted by computational methods. Here, we demonstrate that similar approaches can be used to forecast CRISPR/Cas9 gene editing outcomes in Xenopus tropicalis, Xenopus laevis, and zebrafish. We show that a publicly available neural network previously trained in mouse embryonic stem cell cultures (InDelphi-mESC) is able to accurately predict CRISPR/Cas9 gene editing outcomes in early vertebrate embryos. Our observations can have direct implications for experiment design, allowing the selection of guide RNAs with predicted repair outcome signatures enriched towards frameshift mutations, allowing maximization of CRISPR/Cas9 phenotype penetrance in the F0 generation.Research in the Vleminckx laboratory is supported by the Research Foundation—Flanders (FWO-Vlaanderen) (Grants G0A1515N and G029413N), by the Belgian Science Policy (Interuniversity Attraction Poles—IAP7/07) and by the Concerted Research Actions from Ghent University (BOF15/GOA/011). Further support was obtained by the Hercules Foundation, Flanders (Grant AUGE/11/14) and the Desmoid Tumor Research Foundation and the Desmoid Tumour Foundation Canada. T.N. is funded by “Kom op tegen Kanker” (Stand up to Cancer), the Flemish cancer society and previously held PhD fellowship with VLAIO-HERMES during the course of this work. D.T. and M. C. hold a PhD fellowship from the Research Foundation-Flanders (FWO-Vlaanderen). The Zorn Lab is supported by Funding from NIH National Institute of Child Health and Human Development (NICHD) P01 HD093363. A.W. and A.B. are supported by the Ghent University (Universiteit Gent) Methusalem grant BOFMET2015000401 to Anne De Paepe. The National Xenopus Resource and Horb lab is supported by funding from the National Institutes of Health (P40 OD010997 and R01 HD084409)